Electric-field-induced Mott insulating states in organic field-effect transistors 
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We consider the possibility that the electrons injected into organic field-effect transistors are strongly 
correlated. A single layer of acenes can be modelled by a Hubbard Hamiltonian similar to that 
used for the re-(BEDT-TTF)2X family of organic superconductors. The injected electrons do not 
necessarily undergo a transition to a Mott insulator state as they would in bulk crystals when the 
system is half-filled. We calculate the fillings needed for obtaining insulating states in the framework 
of the slave-boson theory and in the limit of large Hubbard repulsion, U. We also suggest that these 
Mott states are unstable above some critical interlayer coupling or long-range Coulomb interaction. 



I. INTRODUCTION 

Recent advances in producing organic field-effect 
transistors 1 and the growth of single crystals of acenes 2 ' 3 
raises the possibility of using an electric field to inject 
charge into a single layer of a layered molecular crystal 
such as pentacene. Recent theoretical papers have con- 
sidered the possibility of superconductivity due to the 
electron-phonon interaction. 5 ' 4 Superconductivity might 
appear when the filling of the layer at the interface be- 
tween the molecular crystal and the dielectric reaches 
half- filling 4 . In that theory, this would be the conse- 
quence of the existence of a van Hove singularity in the 
density of states and the pairing is described by the 
conventional BCS theory involving the electron-phonon 
interaction. 4 However, in other organic materials, such 
as the k-(BEDT-TTF) 2 X family, the proximity of the 
superconducting phase with various other ground states 
(a Mott insulating state for example 6 ) may suggest a 
more exotic type of pairing, involving the Coulomb in- 
teraction between the electrons. 7 ' 8 It is therefore impor- 
tant to study whether the injected carriers could form a 
strongly correlated state (at half-filling for example). 

The first question we have to study is whether a Mott 
insulating state could be induced. Naively, one might 
expect that if the repulsion between two electrons on the 
same molecule is strong enough an insulating state will 
occur at half filling. The effect of disorder on such an 
insulating state in Cgo crystals doped by charge injection 
was recently considered theoretically. 9 Here, even though 
we consider a clean system, it is not obvious that the 
effects of correlation will appear in a field effect geometry. 

In a field-effect transistor, the charges are confined near 
the interface between the molecular crystal and the di- 
electric, because of the strong electric field used to inject 
the charges. The issue of how far extended from the inter- 
face the electronic wave- functions are is very important. 
As the doping increases, it is indeed not clear whether 
the electrons will remain confined to the first layer or 
because of the Coulomb interaction the electrons will oc- 
cupy other states extended far from the interface. There 
is clearly a competition between the Coulomb repulsion 
that pushes the electrons apart and the electric field that 



forces the electrons to stay close to the interface. Then, 
if the electric field is not strong enough, the Mott tran- 
sition will not occur. In the following, we will adopt the 
slave-boson approach to study the Mott transition in the 
presence of the electric field. Slave-boson theory is the 
simplest theory which generically gives a Mott transition, 
and it will allow us to show that such Mott states could 
occur at non half-integer fillings. 

In section II, we will estimate the Coulomb interaction 
for two electrons on the same acene molecule, such as 
pentacene, tetracene and anthracene since these are im- 
portant examples. We will then give a brief review of the 
physics of a single layer of such molecules. However, the 
presence of a finite electric field in field-effect transistors 
allows the electrons to hop onto the neighbouring layers, 
thus modifying the dynamics. In section III, by applying 
the slave-boson theory, we will show that the Mott state 
could occur at larger fillings, that we calculate in the limit 
of strong interaction. In sections III B and III C, we give 
a discussion of the effect of the interactions we neglected 
first, such as the interlayer hopping, and the long-range 
Coulomb interactions. Those interactions could in fact 
destroy the Mott state, and we will calculate the critical 
strength above which the Mott state is no longer stable. 



II. A STRONGLY CORRELATED MODEL FOR A 
SINGLE LAYER IN ACENE MOLECULAR 
CRYSTALS 

In this section, we consider a model for a layer of acene 
molecules. Band structure calculations have shown that, 
for large electric fields applied to pentacene 4 or Ceo (Ref. 
10) field-effect transistors, the electrons are confined to a 
single layer. 

For acene-based materials, we now argue that a sin- 
gle layer can be described by a Hubbard model on 
an anisotropic triangular lattice. Such a short-range 
Coulomb interaction is more valid in the strong doping 
regime where the screening of the Coulomb interaction by 
the other electrons is more efficient. This is the regime 
we are interested in. The Hamiltonian is written: 
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H = - E *y( c i )CT Cj, CT + h.c.) + m, T ni,i (1) 

i 

where C; CT creates an electron [hole] at site i in the lowest 
unoccupied molecular orbital (LUMO) [highest occupied 
molecular orbital (HOMO)], iy are the tight-binding hop- 
ping integrals between molecular orbitals (t^, td 2 , t a ), as 
shown in Fig. 1, and U is the Hubbard repulsion for two 
electrons on the same molecule that we will estimate be- 
low. The kinetic part of the Hamiltonian can be rewrit- 
ten: 

H o = E e k c L c ^ 

k,<r 

e fe = -2t dl cos(fcidi) - 2t d2 cos(k 2 d 2 ) 

- 2t a cos(fcidi + k 2 d 2 ) (2) 

where ta 12 , t a are the hopping parameters in the d X2 ,a 
directions. For example, t dx = 59 meV, td 2 = 88 meV, 
t a = 47 meV for the LUMO orbital of pentacene. 11 

We note that this band structure is similar to that 
for the k-(BEDT-TTF) 2 X family 8 . A similarity between 
the band structure of a-sexithiophene and this family has 
been pointed out previously. 12 




a 



FIG. 1. Model for a layer of pentacene molecules, accord- 
ing to Ref. 11, 13. The thick lines represent the pentacene 
molecules. t a , ta x and td 2 denotes the different hopping inte- 
grals between neighboring molecules. 
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FIG. 2. Fermi surfaces for the band structure given in equa- 
tion (2), for different band fillings, 1/8, 1/4, 1/2. 



Importantly, the coupling between the layers has been 
found to be negligible, thus showing the strongly two- 
dimensional character of these materials. 11 . There are 
two molecules per unit cell but the band structure re- 
duces to one band due to the equivalence of the crystal 
environment up to a rotation. 

Correlation effects are known to be important for sin- 
gle isolated acene molecules, as has been discussed in 
the quantum chemistry literature 15 ' 16 . For the purpose 
of estimating the Hubbard repulsion, we have calculated 
the Hiickel energies and wave- functions of the electrons 
in benzene, naphthalene, anthracene, tetracene and pen- 
tacene molecules. We have then estimated the cost in 
Coulomb energy of two additional electrons (or holes) on 
the molecule, 

U = J d 3 ri d 3 r 2 ^ (ri)^i (ri) ^ ^ <f* (r^i (ft) (3) 

where <j) a (r) is the lowest unoccupied molecular orbital 
(LUMO) of the molecule. Using the Hiickel approach, we 
found the coefficients q defined by: 

i<m = -4E c ^< ct > ( 4 ) 

* i 

where \i) is the atomic orbital of the zth carbon atom. 
Inserting this into (3) leads us to estimate the following 
expression for U : 

i 

where the expression is written as a function of the lo- 
cal Coulomb interaction on a single carbon atom Uq = 
(i, T;i)J. \e 2 /\ri — r 2 \\i, T; i, I)- We have calculated the 
Hiickel parameters Ci for the four isolated molecules. For 
the benzene molecule, the electron is uniformly delocal- 
ized on the ring. The effective U is therefore Uo/N. For 
the other molecules, there is incomplete derealization 
and £iM 4 > £Jci| 2 = N. Therefore U > U /N. 
The explicit calculation gives U/Uq = 0.14 for naphtha- 
lene (N = 10), 0.12 for anthracene (N = 14), 0.10 for 
tetracene (N = 18) and 0.09 for pentacene (N = 22). 
Assuming a value of Uq = 12 eV for the Coulomb inter- 
action on a carbon atom, 15 one can estimate that U = 1.7 
eV for anthracene, U — 1.2 cV for tetracene and U = 1 
eV for pentacene. Although this calculation neglects the 
screening effects due to the other electrons, it suggests 
that the Coulomb interaction (U ~ 1 eV) may be com- 
parable to the kinetic energy (with a bandwidth given by 
W ~ 0.5 eV 11 ). This is consistent with other similar or- 
ganic materials based, for example, on the BEDT-TTF or 
BETS molecules. 8 Although our estimate of U is rather 
crude (especially because it neglects screening effects), it 
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seems plausible to view a single layer in isolation as a 
strongly correlated system. 

We first review the physics of the one-layer model at 
half-filling 8 and low temperatures before considering the 
effects due to the next layers. At small U <C t, the sin- 
gle layer system, described by the Hamiltonian (1), is 
metallic, because the frustration of the lattice gives an 
imperfect nesting of the Fermi surface (see Fig. 2). At 
large U 3> t, the system is a Mott insulator and has 
a charge gap U. The interaction between the spins de- 
fines a smaller energy scale, <~ t 2 /U. In the intermedi- 
ate regime of U <~ t, the question of the nature of the 
ground state remains to be settled. As previously sug- 
gested on theoretical grounds, 17 a superconducting phase 
may appear near the mctal-Mott insulator boundary. 
This may explain the phase diagram of half-filled organic 
materials. 17,8 As the model (1) is essentially the same as 
the one which describes k-(BEDT-TTF) 2 X, 8 it is tempt- 
ing to suggest that superconductivity may appear in 
acene field-effect transistors, at half-filling. In such an ap- 
proach, the T c is expected to increase with the Coulomb 
interaction. 18 Meanwhile, the larger the molecule, the 
smaller the U is. Therefore, a larger T c is expected for 
smaller molecules such as tetracene and anthracene (if 
screening effects can be ignored). However, the same 
trend is found within the electron-phonon mechanism. 5 
If such strongly correlated superconducting phases ap- 
pear, Mott insulating states should also appear in other 
acenes with a larger U/t. This could be achieved in prin- 
ciple in materials based on a smaller acene molecule (to 
increase U), such as naphthalene. Alternatively, t might 
be decreased by pushing the molecules further apart by 
intercalation or attaching side groups to the molecules. 19 



III. A MODEL FOR LAYERED FIELD-EFFECT 
TRANSISTORS 

We now examine how the single layer approach is mod- 
ified when one takes into account the existence of the next 
layers. We have discussed above the example of acene- 
based materials to give an idea of the order of magnitude 
of the parameters, but the following considerations do 
not actually depend on the details of these materials. 

Our aim is to study the case U 3> t to see whether the 
Mott state we have discussed for a single layer can re- 
ally be induced when the system has many layers. As it 
turns out to be the case, we may then speculate that su- 
perconductivity in field-effect transistors would probably 
be of the same origin as that observed in bulk crystals 
of k-(BEDT-TTF) 2 X. Conversely, the existence of su- 
perconductivity in those materials should encourage the 
search for superconductivity in field-effect transistors. 

According to the band structure calculations, 4,10 the 
electrons are confined to a single layer for the largest 
electric fields applied. Nevertheless, we do not restrict 
ourselves to a single band for two reasons. First, even 



when the charges are confined to a single layer, the next 
layers offer virtual states that make the system different 
from bulk crystals. This is precisely what we study in 
the next paragraph, where the interlayer coupling is set 
to zero, so that there is no mixing of the bands. Second, 
when the electric field is small enough with respect to the 
interlayer coupling t± , the charges are partly delocalizcd 
onto the next layers. So the question is then: can we still 
have a Mott insulating state in the first layer? 

We now take into account the next layers away from 
the interface. The Hamiltonian includes the kinetic en- 
ergy of the electrons within a layer, the strong interaction 
between them, the electric field imposed by the gate (but 
we neglect the modification of the electronic orbital by 
the electric field) and the interlayer hopping: 

U = H e fc c k,i, CT c M,a + U n i,j,l n i,j,l 

+ A ^ i n KhrT - ti^ (4,i J<T c M+i,<r + h- c -) (6) 



where c! . „ creates an electron in a state of momentum 
k, in the layer i and spin state a. is the electronic 
dispersion within a layer, characterized by a bandwidth 
W. A = eEc is the energy spacing between the orbitals 
which belong to different layers, due to the electric field E 
imposed by the gate, c is the distance between the layers. 
We also define 7 = 2A/W . The Hamiltonian is then a 
multi-band Hubbard model with the set of parameters: 
U/W, 7, t±/W and the total filling n. 

dielectric 




first layer 
A second layer 

2A 

FIG. 3. Schematic picture for the first layers near the in- 
terface with the dielectric crystal. A = eEc is the energy 
associated with the electric field. 



A. No interlayer hopping 

Firstly, we will consider the case of decoupled layers, 
t± = 0, which is in fact a good approximation for acene 
materials. The only additional ingredient compared to 
section II is the existence of additional states in the next 
layers at energies [A* — W/2, Ai + W/2] (where A; = iA), 
which possibly overlap with the first band. 
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FIG. 4. Non-interacting band picture with the first two 
bands for a total filling of one electron per molecule, i.e., a 
total of N carriers where N is the number of sites within a 
layer, W the bandwidth and A = eEc the mid-energy of the 
conduction states of the second layer, (a) A > W/2. The 
lower band is half-filled, (b) A < W/2. The lower band is 
no longer half-filled : some electrons occupy the states of the 
second band. 

We now discuss the effects of these additional states 
for both the case of half-filling and arbitrary fillings. 

Half-filling. As long as A > W/2, no additional states 
are filled (Fig. 4(a)) and due to the absence of coupling 
between the bands, the system behaves as a single layer 
in isolation, i.e., there is a metal-insulator transition at 
finite U , as we explained above. The only difference with 
the single layer problem is the existence of additional 
excited states, which will affect the transport and mag- 
netic properties of the system. Thus, depending on the 
strength of U/A, the system is either a Mott insulator 
(U/A < 1, the charge gap is roughly given by U) or 
a charge-transfer insulator (U/A > 1), for which the 
excited states consist of a charge transfer to the next 
layer (the charge gap is given by A). Thus, the gap can 
be much smaller than U if the electric field is not large 
enough. The magnetic properties of the layer are given 
by superexchange processes. In the case of the Mott in- 
sulator, the effective antiferromagnetic exchange is given 
as usual by J = At 2 /U. In the case of the charge-transfer 
insulator, a given electron would not necessarily hop vir- 
tually to a neighboring site in the same layer because it 
costs the energy U but would rather hop onto the next 
layer which costs a smaller energy A = eEc. Then, in this 
case, the exchange J between two spins would explicitly 
depend upon the strength of the electric field. Note, how- 
ever, that if we change E, we would not only change J, 
but also the total filling of the system. Above half-filling, 
some carriers will occupy states in the next layer. As 
long as the layers are weakly coupled, these extra carriers 
will not change the magnetic properties of the first layer. 
Therefore, the net effect of the electric field will be to de- 
crease the interaction between the spins. This effect can 
not be detected through the change of the critical tem- 
perature of the antiferromagnetic transition because such 
transition is prevented by the Mcrmin- Wagner theorem 
in a strictly two-dimensional system. However, measure- 



ment of the magnetic susceptibility would provide infor- 
mation about spin-spin interaction and how they change 
with the electric field. When A < W/2 (see Fig. 4(b)), 
however, the second layer starts to be filled before the 
first one gets to half-filling (the chemical potential is neg- 
ative). The entire system is less than half-filled and no 
metal-insulator transition is then expected. The critical 
value A = W/2 is very schematic here because t± = 0. 
In the section IIIB, we will see that for finite t±, the 
critical value of A actually increases. 

Arbitrary fillings. In the case A < W/2, some states 
in the second band have a negative energy and need to 
be filled before the first band becomes half-filled (see 
Fig. 4(b)), thus possibly allowing for a localization of 
the charges in the first layer if U is large enough. To as- 
sure that the first band is half-filled, we fix the chemical 
potential at zero energy, /i = 0. In the next paragraph, 
we calculate, in the framework of the slave-boson theory, 
the total number of carriers needed to have /j, = 0. 

In the slave-boson mean-field theory, 20 the effect of 
the Coulomb interaction is to renormalize the hoppings 
by quantities denoted by Zj (for the layer j). Zj is the 
mean field value of the operator zij which is given by 
a combination of auxiliary boson operators, e\jPi,j,a + 
p\ j _ cr dij , that tell if the site is empty (operator 

eij), single occupied (pi,j, a ) or double occupied (dij). 14 
The many-body Hamiltonian (6) for the layer j reduces 
in the mean-field approximation to free fcrmions with 
renormalized bandwidths: 

il.cr i,<7 

(7) 

where dj is the mean-field value of the double occupancy 
operator dij. In the following, we will restrict ourselves 
to U —* +oo, so that dj — (the double occupancies are 
forbidden). The layers are all decoupled and the total 
Hamiltonian is the sum over all the layers. However, 
the quantities z 2 , which determines the bandwidths, arc 
coupled together because they depend on the fillings of 
the band j, which is determined by the chemical potential 
of the entire system. 

For free electrons (U = 0, z| = 1), the number of states 
below zero energy in the second layer is denoted by (1 — 
j)N (7 < 1 and N is the number of sites). When we take 
into account the Coulomb interaction within the slave- 
boson picture, the bandwidths are reduced (W — ► Wz\). 
The number of states given above becomes (1 — 7/z 2 )7V, 
that is smaller than (1 — j)N and depends upon z\ that 
we have to determine. Therefore, we note that the total 
filling needed to have \x = in the interacting system 
must be smaller than that for free electrons. We now 
calculate the filling explicitly. 

We first calculate the value of z\ for 7 < 1 when the 
carriers of the first layer become localized. For that, we 
consider just the first two layers (j = 0, 1) and take a 
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constant density of states in the layers to simplify the 
calculation. At /j = 0, the fillings nj are given by: 



1 (n 

2 + ~ 



i)4 



■7/2 



z 2 + zl 



n 2 = n — ni 



(8) 
(9) 



where n is the total filling of the system, Wz\ and Wz\ 
are the bandwidths of the two bands. In the limit of large 
U, (i.e., U > W, A), z\ = (1 - 2m)/(l - m). Charge 
localization in the first layer is characterized by z\ = 0, 
or m = 1/2. First, the equation (8) gives the bandwidth 
of the second band z\ — 7/[2(n — 1)]. On the other 
hand, in the limit of large U, the bandwidth is given by 
z\ = (1 — 2ri2)/(l — ri2), where n 2 = n — 1/2. Equating 
the two expressions, we find the total filling n, n c (j), for 
which m = 1/2 and the bandwidth z%: 
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n c {j) = 1 - 



87- 
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7) 



V 1 + 8/7-1 



= 1 - /(7) 



(10) 
(11) 



Note that this solution is valid as long as the third band is 
not involved, that is to say, if we define the energy of the 
middle of the third band by A 2 = a 2 A, A 2 - W/2 > 0, 
or in terms of 7, 7 > l/a 2 . 

Therefore, in the regime l/a 2 < 7 < 1, the total filling 
needed to get the first band half- filled is given by eq. (11). 



Above this filling, 



0, so that we have one localized 



electron on each site of the first layer. The additional 
(2n c (7) — 1)N electrons fill the band of the second layer 
which bandwidth is W 'z\ (given by eq. (10), see Fig. 5). 



n<n (y) (a) E 



Wz 2 2 A 
-tWzfO 




n>n c (y) 



(b) 



Wz 2 2 
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FIG. 5. Interacting band picture for the first two 
bands within the slave-boson mean-field theory (here 
2 A < W < U). (a) The first two bands are partially filled 
with electrons and the bandwidths are renormalized by the 
interaction, (b) As soon as n > n c (7), the lower band is 
half-filled, so that z\ = 0. Then, there is one electron lo- 
calized on each site of the first layer. The second band is 
partially filled. 

More generally, in order to consider the regime with 
7 < l/ajvt, we have to consider Nf, additional bands. 
Generalization to many bands with energies Aj = o^A 
(ao =0, ot\ = 1, a.2---) and bandwidth W is straightfor- 
ward. Performing the same steps as above, in the regime 
l/a^+i < 7 < 1/ajVb, the filling at which the Nb first 
bands arc half-filled is given by: 



N b + l 



(12) 



P =i 



Let us recall that Nf, is the number of bands partially 
filled. For a given system with a fixed 7, N}, bands are in- 
volved when the bottom of the N^ h band crosses the zero 
energy. It occurs when 7 gets smaller than 1/a^.The 
shape of this curve depends on the positions of the bands, 
Ap = apA. 

Furthermore, the generalization to higher fillings is 
also very similar. If we consider such a filling that P 
bands are completely half-filled and Nb bands are par- 
tially filled. In the regime where l/(ap+7v b — ctp—i) < 
7 < l/(ap + N b+ i — ap-i), the critical n c to get the P th 
band half-filled is given by: 



rich) = 



N b 



N b +P-1 
3=P 



Note that this equation reduces to (12) for P ■■ 
amples of the first two curves (P = 1 and P 
given in figure 6. 



(13) 

1. Ex- 

2) are 



O 

0.5 




the first two layers 

are half-filled (P=2) 



the first layer is half-filled 

(P=D 



u =< 



0.5 



7=2A/W 



1.5 



FIG. 6. Phase diagram (filling, n, versus energy spacing be- 
tween the bands due to the electric field, 7) from slave-boson 
theory in the limit U = 00 and t± = 0. Solid lines show the 
critical fillings above which the first (resp. second) band is 
half-filled, i.e. the electrons are localized on each site of the 
first (resp. second) layer in the limit of large U. The dashed 
line shows the critical filling above which the second band 
starts to be filled. 

In conclusion, when 7 > 1, only one band is involved. 
When the band gets to half-filling (n = 1/2), it under- 
goes a transition to a Mott state at large U in which all 
the electrons are localized on the first layer. Doping the 
system further leads to filling of the states of the second 
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layer, eventually up to half-filling when the total density 
is n = 1, then localizing the electrons in the second layer, 
and so on (see Fig. 6). However, as soon as 7 < 1 (the 
second band starts to be involved above some filling, see 
dashed line in Fig. 6), the Mott state needs a larger fill- 
ing to occur. This filling, n c ("f), is given by eq. (12). In 
that case, the first band is half-filled, and there is one 
electron localized on each site. As the total number of 
electrons is larger than the number of sites of the first 
layer, some of them occupy the states of the next layers. 
The system is metallic, but with a reduced number of 
carriers. 

The condition on 7 is nothing but a condition on the 
electric field since 7 = 2eEc/W. The border 7 = 1 
defines a critical field E c . For instance, in the acenes, 
W ~ 0.2 eV and c - 15A gives E c ~ 10 s Vm" 1 , which 
is close to the dielectric breakdown field. 21 

However, this calculation neglects the interlayer hop- 
ping. It is known that such couplings arc quite small 
in layered molecular crystals, in particular in those dis- 
cussed above, 11 ' 13 ' '' 19 but it is not obvious whether an in- 
finitesimal coupling would destroy the charge localization 
in either layers. In the next section, we study whether 
such a state survives when one takes into account such 
interactions, such as the interlayer coupling or also the 
Coulomb interaction between nearest sites. 



B. Effect of the interlayer hopping 

The Mott state in a single layer of a field-effect tran- 
sistor is not as stable as in a bulk crystal. The elec- 
trons have, indeed, other available conduction states in 
the next layers away from the surface. This opens a way 
to avoid the strong on-site Coulomb interaction, though 
it costs the energy associated with the electric field, eEc. 
Some interactions will actually act against the charge 
localization: the coupling between the layers favors ex- 
tended states away from the surface. We just consider the 
hopping between the two layers closest to the surface: 

W± = -U^(4 ;liCT c k , 2 , a +h.c.) (14) 

k,cr 

When 7 > 1, the two bands are separated by a gap. 
Using the same Hamiltonian, but at half-filling, Blawid 
et al. showed that, within the slave-boson theory, there is 
a Mott-metal transition when t± > W[y(-y-l)/2\ 1 / 2 /2 = 
t± c for U = +00. 22 Note that, when 7 — > 1, the critical 
t± vanishes. This is consistent with the disappearance of 
the Mott state at half-filling that we found when 7 < 1. 

On the other hand, in section III A, we found that a 
total filling larger than n c (^) (> 1/2) was needed, when 
7 < 1, to have the first layer half-filled, thus allowing for 
a Mott state in the first layer. We now wonder if this 
state is still stable above n c (^) when t± 7^ 0. If we think 
of the result given above, it seems that we need a real gap 



between the two bands (7 > 1). With the Coulomb in- 
teraction switched on, the bandwidth of the second band 
is reduced and a gap opens when A - Wz\j1 > 0. In 
the limit of large U, z\ is given by (1 — 2n 2 )/(l — n 2 ), 
and n 2 = n — 1/2 as long as any other band can be 
safely ignored. Therefore, the gap opens when the den- 
sity n satisfies z\ = 7 = (2 — 2n)/(3/2 — n). That 
gives n = (2 - 37/2)/(2 - 7) = n' c (j) > n c ( 7 ). So we 
could think that the Mott state is unstable in the regime 
n c (7) < n < n' c (j), but becomes stable as soon as a gap 
opens, i.e. when n > n' c (j). 

In fact, the case considered here is different from that 
of Blawid et al. because the second band is partially 
filled. We now show that, as soon as n > n c (j), the 
Mott state is stable below a critical value for t±. This 
is because we do not need a gap between the two bands, 
but a gap between the energy of the localized states and 
the energy of the available states in the second band. As 
soon as \i > 0, the second band is partially filled and the 
gap opens. 

We consider two coupled bands with dispersions given 
by efe and (?& = A + in the non-interacting pic- 
ture) . We neglect the other higher-energy bands in order 
to simplify the expressions. In the slave-boson picture, 
the effect of the interaction is just to rcnormalizc the 
bandwidths by quantities denoted by z\ and z\. 

H o = z hkct^ a c Khr7 + ( A + z hk) cL, ff c k,2,<T (15) 

k.cr k,/ ? ' 

It also renormalizes the coupling between the two bands 
such as: 

H± = -t±Z!Z 2 ^2 (Ck,i,a c k,2,CT + h.c) (16) 

k,i,(T 

We again restrict ourselves to the limit of infinite U, 
where the zf are function of the fillings rij: 



»i=2^E< C U^ fc ,i, < r> (18) 

where the m are calculated as function of the parame- 
ters of the model. Near the metal-insulator boundary for 
the first layer, (but in the metallic region), m will be 
expanded as a function of z\ — > 0: 

ni = \ - Pzl (19) 

where we have introduced a parameter fj that has to be 
determined. To satisfy equations (17) and (19) at the 
same time, in the metallic regime, we need: 

4/3 = 1 (20) 
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This gives a condition on the parameters of the model 
at the Mott-metal boundary. So all we have to do is 
to calculate this parameter (3. We decompose c\ 1 a as 
function of the eigen-operators: 



J 

-k,l,a 



(21) 



where u k and v k are given by the diagonalisation of the 
Hamiltonian. We are interested in the limit t± <~ z\ — > 0. 
Then, these parameters are expanded as: 



\Uk\ 



1 



ti 



(efe - hf 



(22) 
(23) 



and m can be rewritten: 



fe fe 



i. A* = 



(24) 



We restrict ourselves to 7 > 1, as the Mott state ap- 
pears to be unstable when 7 = 1. Then only the lower 
band is filled, i.e. ((3l a Pka) — and m is given by: 



= ---E- 



- \2 



(25) 



where the sum runs over all the k. At the first order in 
zf, we can replace e k by in the expression above, e k has 
a bandwidth which depends on z\. As the second band 
/3 is empty, ni is of the order 0(tjj. Then, using (17), 
z\ is just 1 — 0(t\). We can replace i k — e k by A + e fe : 



ni 



-E 



t 2 



2N^(A + e k ) 2 



(26) 



which can be expressed as a function of an integral: 



1 1 

" 1 = 2 _ 2 



W/2 



dujpo((jj) 



(A + u;)< 



(27) 



with po = 2/W, the density of states per site that we have 
taken constant. The upper limit is 0, because only half 
of the first band is filled. Now, we make the substitution 
t± — > t±Z\: 



therefore, the critical t±, given by eq. (20), is: 



t ±c = VA(A - W/2)/2 



(29) 



which is the expression found previously. 22 Along the line 
7 > 1, n = 1/2, the Mott state is stable below a critical 
t± which is given above. Note that it vanishes when the 
gap between the two bands vanishes, 7 — > 1 (A — > W/2). 



2. p > 

We have to fill the second band (5 with electrons, m 
has an additional term, and is given by: 



2 2N 4^ 



1 



2]y "^(fikaPka) 



(30) 



where {0l a l3ka) = 0(i k - P + We have now to 

calculate the third term. The second one is, indeed, iden- 
tical, provided that W — > Wzf in order to take into ac- 
count the renormalizcd bandwidth of the second band 
due to its filling. We first calculate the integral and then 
make the substitution W 



Wzl 



J -W/2 



ti 



dojp (uj)- 
W/2 (A + w) 2 



t\ p - A + W/2 



W p(A - W/2) 
Therefore, after the substitution, we find: 



(31) 



n l = 77 



ti 



■ + 



t\ p-A + Wzl/2 
2A(A - Wzl/ 2) ' Wzl p(A - Wz 2 /2) 

The next question is how p is related to n. 



- m = - / <kjp(us) 

1 J -W/2 



(32) 



(33) 



where p(u) is the density of states of the second band, 
which is affected by the coupling t± . At the leading order, 
however, we can replace p(uj) by pq{oj) and n\ by 1/2. 
That gives the chemical potential: 



p = A + (n - \)Wz\ + 0{t\) 



(34) 



The problem is that z\ is also function of the filling n 2 . 
But here again, because we are looking for the leading 
order in z 2 , we can replace n 2 by n — 1/2. Then, 



t 2 z 2 



2 2A(A - W/2) 



(28) 



7 2 _ 
z 2 — 



1 - 2n 9 



2-2n 



1 - n 2 3/2 - n 



(35) 
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Putting all the terms together and noticing that t± scales 
as Z1Z2, n\ is then given by: 



t±zi 
W 



3/2- 



2n-l 



7- 



2-2w 
3/2-n 



7 + 2(n-l) 



2-2w 
3/2-n 



and the equation for the stability border, 4/3 c = 1, is 
written as: 



t±c 

W 



4— An 
3/2-n 



In- 1 



7' 



2-2n 
3/2-n 



7 + 2(n- 1) 



3/2- 



-1/2 



(36) 

The resultant phase diagram is given in Fig. 7(a) for a 
particular value of 7 < 1. Below the critical t± Cl the first 
layer is insulating. The t± c vanishes at a particular den- 
sity. As shown in Fig. 7(a), it vanishes with /j, , which 
corresponds to the density n = n c (S) that we calculated 
in the previous section (for the present problem of two 
bands, see eq. (11) or eq. (12) or (13) when more bands 
are involved). In other words, t± c vanishes when there 
is no gap between the energy of the localized states and 
the energy of the available states in the second band. In 
conclusion, as soon as n > n c (5), there is one electron 
localized on each site of the first layer for t < t± c . Re- 
member, however, that we are considering the U = 00 
limit. 

Note that, when n = 1/2, the second term in the brack- 
ets vanishes, leading to the expression (29) which is valid 
for A > W/2 (7 > 1). 

In Fig. 7(b), the critical t± c is given for several values 
of 7. When 7 > 1 (see the case 7 = 1.2), there is a 
critical t± even at n = 1/2. When n increases above 1/2, 
the bandwidth of the second band decreases according to 
eq. (35) and the gap increases. Consequently, the mixing 
of the two bands decreases and we need a larger t± c to 
destroy the Mott state. At n = 1, it diverges. In the 
model we are considering, we have only two bands. When 
n = 1, the entire system is half-filled and is therefore a 
Mott insulator. The tj_ c to destroy it is then infinite 
because we are in the limit of infinite U. 





(a) 




0.4 


/ 7 = 0.8 : 




Metal 


/ Insuidtor 


\ 
■+-> 

0.2 






0.1 




/m/"w 









1(b) 


1 1 1 r , 1 1 1 r | y , . i , i , y i 


- 7/1.2 / / : 




y 1 / 




/ 7 7°- 8 ■ 




j 1 n o(7) ; 



0.4 0.6 



0.45 0.5 0.55 0.6 0.65 0.7 



FIG. 7. (a) Example of the phase diagram when t± / 0, 
for 7 = 0.8. Charge localization on each site of the first layer 
occurs as soon as the chemical potential n > 0, i.e. as soon as 
there is a gap between the energy of the localized states and 
the energy of the available states in the second band, (b) The 
critical line for different values of 7. In every case, U = 00. 



C. Inter-site Coulomb Interaction 

Another example of an interaction that can destroy the 
Mott state is given by the Coulomb interaction between 
the nearest neighbor sites: 



(37) 



<i,i> 



This term is diagonal in the basis of localized electrons. 
In the limit of large U and at half-filling, the energy per 
electron of the Mott insulating state is Eo = zV/2, where 
z is the coordination number of the lattice. Therefore, we 
expect that if E becomes larger than the energy spacing 
between bands, A, it will be energetically favorable to 
transfer carriers from the Mott insulating layer to the 
next band. The system will then be metallic. Therefore, 
in the limit of infinite U, we need V < 2A/z to get a 
Mott state. 



IV. CONCLUSION 

We have given an estimate of the on-site Coulomb in- 
teraction for the layered acene molecular crystals. It 
turns out that the interaction (U ~ 1 eV) may be compa- 
rable to the band width (W ~ 0.5 eV). However, whether 
the system is in the metallic or insulating regime, has to 
be determined experimentally. A prediction of the na- 
ture of the real ground state would indeed need a more 
accurate estimation of U/W. This is in fact a common 
feature of many organic materials 8 . 

In bulk crystals, a large repulsive interaction leads to a 
Mott insulating state when there is one electron per site. 
In such a case, the hopping of the electrons to the other 
sites costs the large energy U. In field-effect transistors, 
however, when the first molecular layer at the interface 
reaches half-filling, the electrons still have many empty 
states available on the next layers away from the surface, 
though these states are at higher energies depending on 
the strength of the electric field. When the electric field is 
smaller than a critical value, E c — W/(2ec), the carriers 
occupy the conduction states of the next layer and the 
first layer is doped with holes. In this case, the Mott 
state is destroyed. However, above the critical electric 
field, the Mott state is stable at half-filling, at least below 
a critical value for the interlayer coupling, t± c (given by 
eq. (29)). 

On the other hand, we have shown that it is possible to 
restore the Mott state below E c provided that the system 
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is further doped. This is in order to fill the first band up 
to half-filling, to compensate for the loss of carriers which 
go onto the next layers. At n = n c (7), the chemical po- 
tential vanishes and allows for a Mott state in the first 
layer. This is what we have found in section III A, ne- 
glecting the coupling between the layers. This coupling, 
although small in the layered molecular crystals, is still 
present. In section IIIB, we have calculated the criti- 
cal value, t± c , below which the Mott state is stable. It 
turns out that as soon as n reaches n c (-f), t± c becomes 
non-zero. 

In reality, in a field-effect transistor, the doping n and 
the gate voltage (or equivalently the electric field E) are 
not independent quantities. On Fig 7(b), if we take a 
particular material with a fixed t± and W, when the 
electric field increases, the metal-Mott boundary will be 
crossed at some stage (provided that U is large enough) . 

Finally, we speculate about the possibility of super- 
conductivity in organic field-effect transistors. Since a 
Mott state can be induced in the first layer if the inter- 
action is large enough, it is possible that superconductiv- 
ity can occur due to the strong correlation between the 
electrons (though not strong enough to create a Mott 
insulating state). Is it possible to find a superconduct- 
ing phase near the metal-insulator boundaries as in n- 
(BEDT-TTF) 2 X? Then, if the superconductivity is in- 
duced in the first layer, does it effect the next layers by 
the proximity effect? Experimentally, a strong electric 
field (E ~ W I (ce)), is needed, either to reach a sufficient 
doping or to get a sufficiently large energy separation be- 
tween the bands. It is not obvious whether the regime in 
question can be experimentally reached or not because 
of dielectric breakdown; but it turns out that layered or- 
ganic molecular crystals are good candidates to observe 
such effects because of their relatively small bandwidth 
and weak interlayer coupling. 
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